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Abstract 

We advocate the replacement of standard a s (/x)-based QCD perturbation 
theory, in which the coupling and truncated perturbative predictions are 
dependent on the chosen renormalisation scheme, by a A-based approach in 
which QCD observables are directly related to the dimensional transmutation 
parameter, A, of the theory. The shortcomings of the standard approach are 
emphasised by formulating it in the A-based language. We show how the A- 
based approach can be extended to include massive quarks. We also consider 
the use of A-based perturbation theory in the case when a resummation of 
large infrared logarithms is required, performing a phenomenological analysis 
for the e + e~ event shape observables thrust and heavy jet mass. We show 
by fitting to data that smaller power corrections are required than in the 
standard approach. 
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1 Introduction 



The process of removing ultraviolet divergences in defining the renormalised 
coupling of quantum field theories is not uniquely specified since the finite 
part may be adjusted at will. This renormalisation scheme (RS) dependence 
means that renormalisation group improved (RG-improved) fixed-order per- 
turbative predictions depend on the arbitrary renormalisation convention 
used. The small value of the coupling in Quantum Electrodynamics means 
that this is not a practical problem. As is well-known, perturbative calcu- 
lations with the conventionally used on-shell renormalisation scheme are in 
exceptionally good agreement with experimental measurements of, for in- 
stance, the anomalous magnetic moment of the electron. The situation in 
Quantum Chromodynamics is not so favourable and the scheme dependence 
problem is a source of serious ambiguity in making precise perturbative QCD 
predictions, for a review see pQ. We shall refer to the standard approach as 
"a s (/i)-based" perturbation theory (a s (/z)BPT). The renormalisation scale is 
conventionally chosen to be fi = xQ, where Q is the physical energy scale of 
the process, for instance the centre of mass energy a/s in e + e~ annihilation. 
The modified minimal subtraction (MS) procedure 2\ is most commonly 
used where the ln(47i~) —Je which always arises using dimensional regulariza- 
tion is subtracted off. Both the choice of the dimensionless number x in the 
scale choice and the subtraction procedure are of course completely arbitrary. 
Worryingly for many QCD observables calculated perturbatively at next-to- 
leading order (NLO) in perturbation theory, the unphysical x-dependence of 
the perturbative prediction is quite strong. The conventional attitude to this 
is to vary x around the so-called "physical scale" /i = Q , i.e. x — 1, be- 
tween say x = 1/2 and x = 2, the resulting a s (xQ) values are then evolved to 
Q = Mz using the QCD beta-function equation truncated at the appropriate 
order, and a central value of a s (Mz) corresponding to x — 1 together with 
a "theoretical error bar" corresponding to the variation of x over the cho- 
sen range, are obtained. The difficulty with this approach is that the range 
of scales to be considered is completely ad hoc. There is nothing sacrosanct 
about the choice fi = Q. Indeed the meaning of fi depends on the subtraction 
procedure employed. For instance the choice fi = Q with MS subtraction 
gives exactly the same result as fi = 0.38Q with MS subtraction where the 
In 47r — 'Je is n ot subtracted off. Clearly the choice of scale cannot be phys- 
ically motivated. Unfortunately there are many arguments in the literature 
which attempt to do just this. A recent example is Higgs production by bot- 
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torn quark fusion where a kinematical argument is used to infer the correct 
choice of the factorization scale For truncated fixed-order perturbative 
predictions it was shown in jl] how the RS-dependence of the renormalised 
coupling and the perturbative cofficients can be parametrised by the com- 
bination [i/Ars, and by the non-universal coefficients of the beta-function. 
Here Ars corresponds to the dimensional transmutation parameter with a 
particular choice of subtraction procedure for removing UV divergences. The 
beta-function is conventionally taken to be truncated at the same order as 
the perturbation series. Whilst the truncated results depend on RS, the 
formal all-orders sum of the perturbation series does not. It was therefore 
advocated in Ref . jlj that the RS in each order should be chosen according to 
the Principle of Minimal Sensitivity (PMS), that is so that the approximant 
is stationary with respect to variation of the scheme parmeters. Whilst this 
idea can be motivated by a number of persuasive toy examples one cannot 
make any firm statements about the goodness of approximation of the PMS 
fixed-order results. One is making the best of the RS-dependence problem 
with which standard fixed-order RG- improved a s (/i)BPT is afflicted. 

In this paper we wish to propose an alternative A-based formulation of 
perturbative QCD (ABPT). The idea is extremely simple, and follows from 
dimensional analysis. It is very closely related to the discussion of [Sj, and 
is essentially the method of effective charges introduced by Grunberg [UJ [7] , 
although as we shall explain, ABPT can, at least in principle, be defined 
even if the QCD observable is not an effective charge. For a dimensionless 
QCD observable 7Z(Q) depending on the energy scale Q, and with massive 
quarks with pole masses {mj} one will have on general dimensional analysis 
grounds, 



Here A is a dimensionful constant related to the dimensional transmutation 
parameter of the theory. We can invert $ to obtain, 



and so the A parameter can be directly related to the measured observable. 
The key point to realise is that the dimensional transmutation parameter 
which is associated directly with broken scale invariance is the fundamental 




(1) 
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parameter. The renormalisation prescription used to remove UV divergences 
will determine the value of A. However A's corresponding to different pre- 
scriptions can be exactly related given a one-loop perturbative calculation 
in the two schemes jH]. Since MS subtraction is most commonly employed 
Ajfg can serve as the fundamental QCD parameter. Clearly the function 
$(1Z(Q), {rrij}/Q) which determines Ajjg does not involve the RS parame- 
ters on which the standard a s (/z)-based fixed-order results depend, and as 
reconstructed from the perturbative information it will automatically involve 
RS-invariant combinations of perturbative coefficients and beta-function co- 
efficients, and will be completely independent of the renormalisation scale p. 
The explicit form of the function $ is easily obtained in the case that 7Z{Q) 
is simply related to an effective charge jUj- An effective charge corresponds 
to an observable having a perturbative expansion, 

TZ{Q) = a{p) + Y j r n a{p) n+ \ (3) 

n>0 

where a(p) = a s (p)/ir is the RG-improved coupling. The coupling a(p) will 
satisfy the beta-function equation, 



da 

— — — p[a) = —ba' 
dm// 



l + ca + ^c n a n ) , (4) 

V n>l J 



where b = (33 — 2Nf)/6 and c = (153 — 19Nf)/12b are the universal beta- 
function coefficients for Nf active quark flavours. If we assume massless 
quarks, so that rrii = 0, then dimensional analysis implies that, 

dH(Q) = p(K(Q) 
dQ Q U 

where p(1Z(Q)) is a dimensionless function of 71, which can be recast as, 

d7l 



dlnQ 



p(K(Q)) . (6) 



The function p(TZ(Q) may be obtained perturbatively by the following al- 
gebraic steps. Starting from the perturbation series in Eq.(3) one chooses 
/i = Q. Differentiating with respect to \nQ term-by-term and using the 
beta-function equation of Eq.(4) one then obtains dlZ/dhiQ as a power series 
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in a(Q). Finally, one can invert the perturbation series in Eq.(3) to obtain 
a(lZ) as a power series in 1Z. In this way one finds, 



p(K) = -blZ 2 ^1 + cJZ + J^PnK"^ 



(7) 



where the p n are RS-invariant and Q-independent combinations of the per- 
turbative coefficients and beta-function coefficients 

p 2 = c 2 + r 2 - cr x - r\ 

P3 = c 3 + 2r 3 - 4r x r 2 - 2r 1 p 2 - cr\ + 2r\ . (8) 



As we shall review in Section 2 one can then integrate Eq.(6), and the A 
parameter arises as a constant of integration. One finally obtains the relation 

An = QF{n{Q))G{1l{Q)) . (9) 

The integration contstant A^ is exactly related to A-j^g (the tilde implying 
the A definition favoured in Ref.[3]) with 

A^ = e^Ajrs , (10) 



Here r denotes the NLO perturbative coefficient r\ of Eq.(3) in the MS 
scheme with renormalisation scale p = Q. J- (TV) is the universal function, 

F(K) = e- 1/m (l + l/cH) c/b , (11) 
and the function Q(7V) depends on the observable and has the form, 

oo 

g(K) = l + J29nK n . (12) 

n=l 

The g n coefficients are RS-invariant and Q-independent combinations of the 
p n in Eq.(8). The first few being, 

P2 

91 = -J 

P3 . P\ P2 n „x 

92 = -T + W~ C T- (13) 
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Eqs. (9), (10) then enable direct extraction of Ajjg from the measured data 
for 7Z. Such direct extractions of Ajjg were carried out for data on e + e~ 
annihilation event shape observables in |5]. The approach is very closely 
related to the method of effective charges of Grunberg [7], and also to the 
Renormalisation Scheme Invariant Perturbation Theory (RESIPE) proposal 
of ^Uj 5 and the Complete Renormalisation Group Improvement (CORGI) 
approach of ^T]. The purpose of introducing the ABPT terminology is to 
emphasise the fundamental importance of the dimensional transmutation 
parameter A, and the dimensional analysis from which it follows. It should 
be noticed that dimensional analysis implies that an equation of the form of 
Eq.(6) holds even if 7Z(Q) is not an effective charge. However, it is not then 
easy to explicitly integrate Eq.(6). Hence, for instance, it is not known how to 
arrive at a result of the form of Eq.(9) relating 7Z and A if 1Z is an observable 
involving the convolution of structure functions and hard scattering partonic 
cross-sections. As we shall review in Section 2 structure function moments 
can, however, be directly related to effective charges. 

The plan of this paper is to begin in Sec. 2 by reviewing ABPT for effec- 
tive charges, as discussed in Ref. In Sec. 3 we show how the standard 
fixed-order a; s (/z)BPT can be written in the ABPT language, this in turn will 
be used to emphasise that the standard approach is very likely to provide 
an erroneous determination of A. In Sec. 4 we show how to generalise ABPT 
to the case of massive quarks, in which case Eq.(6) is no longer a separa- 
ble differential equation. Nevertheless we show that for mass-independent 
renormalisation schemes such as MS Eq.(9) still holds relating 1Z and A, but 
where now the pi invariants of Eqs. (8) are Q-dependent. In Sec. 5 we then 
show how to generalise ABPT to include a resummation of large infrared log- 
arithms in the e + e~ annihilation event shape observables thrust and heavy 
jet mass. We find that significantly smaller power corrections are needed 
to fit the data than with the standard resummations based on the physi- 
cal choice of renormalisation scale. This extends the recent analysis of the 
DELPHI collaboration J2| where ABPT was applied to the analysis of event 
shape means, and where it was found to be possible to obtain consistent 
fits to data without any power corrections. We consider the distributions 
themselves. Finally, Sec. 6 contains a Discussion and Conclusions. 
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2 ABPT for effective charges 



The dimensional analysis result of Eq. (6) is in the form of a separable differ- 
ential equation and may be straightforwardly integrated to obtain 



A 



n 



K(Q) 



dx 
p(x) 



+ K . 



(14) 



The constant of integration which arises is necessarily infinite to ensure that 
asymptotic freedom, TZ(oo) = 0, holds. It can be split into lnA-^ + k,, where 
is a finite, dimensionful and observable-dependent scale, and k is a uni- 
versal infinite constant. The integral in Eq.(14) diverges at x = since from 
Eq.(7) dx/p(x) = l/[x 2 (l + cx + . . .)]. It is then clear that to cancel this 
infinity the infinite part of the integration constant k must be of the form, 



dx 



+ cx + A(x)) ' 



(15) 



where A(x)/x 2 is required to be finite at x — 0, but is otherwise arbitrary, 
and C is an arbitrary upper limit of integration. The choices A(x) =0 and 
C = oo are convenient, but different choices can be completely absorbed by 
a suitable redefinition of A^. Defining k in this way we can regroup Eq.(14) 
in the form 



b In 



Q_ 

A. 



dx 



n(Q) 



[1 + cx) 



X/ 



+ 



rn(Q) 

/ dx 
Jo 



+ 



p(x) x 2 (l + cx) 



(16) 



Denoting the integrals on the r.h.s. of Eq.(16) by F(7Z) and G(1Z) respec- 
tively one finds 



Fiji) = - + cln 



cTZ 



1 + cTZ 



(17) 



Exponentiating yields Eq.(9), with !F{TZ) and Q{TZ) as defined in Eqs.(ll), 
(12). The coefficients gu in Eq.(13) follow on expanding exp[— G(7Z)/b] as 
a power series in TZ. We finally need to relate A^ which is observable- 
dependent, to the universal A^-. This is easily done since on rearranging 
Eq.(9) and taking the limit as Q — > oo one finds 



A^ = lim Qexp{-F{TZ{Q))/b) , 



(18) 
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where the fact that limg^oo G(1Z(Q)) = G(0) = has been used. If Ajjg- is 
obtained by integrating up the MS beta-function equation for a(Q), the MS 
coupling with \i = Q, 

MQ) - BiM , (19) 



cflnQ 

with the same conventions for k, then one similarly finds 

~Ajjs = lim Qexp(-F(a(Q))/b) . (20) 

Noting that 

F(ll)~F(a(Q))-r + ... , (21) 

where the ellipsis denotes terms which vanish as Q —>■ oo, one can then read 
off the exact relation 

A n = e r ' h k- m . (22) 

The tilde over A reflects the fact that the infinite constant k has been de- 
fined with the convention adopted in [I], whereas the commonly assumed 
convention [T3| corresponds to k — > k — cln(b/2c), one then has 

2c^ c/b 



^Ms=yj) Aw- (23) 
We can then arrive at the desired relation between Ajjg and 7Z, 

Aw = QK^T{n{Q))g{n{Q)) , (24) 

where we have defined the observable- dependent normalisation constant 

= e- r ' h {2c/b) c/b . (25) 

Notice that JC^ S is the only part of the expression dependent on the subtrac- 
tion convention used to remove ultraviolet divergences. Changes in subtrac- 
tion convention are trivial in that JC-ji is simply scaled by a constant precisely 
calculable given that the NLO coefficients r are known in the two RS's [S]. 
We have for instance 

'MS _ flE- ln(47T) \ jig 



K%> = exp ^ p^ j K%> . (26) 

Without a NLO calculation the normalisation factor JC-ji is unknown and 
so tree-level calculations by themselves cannot be used to determine Aj^, 
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equivalently in the standard a s (/i)BPT approach tree-level calculations have 
a monotonic \i dependence and are hopelessly ambiguous. Given a NLO 
calculation K-ji is known, but the coefficients gt in Eq.(13) are unknown 
without a N* +1 LO calculation, so that the state of our knowledge is G(7Z) = 1, 
and the deviation from unity is simply unknown without a NNLO calculation, 
which would give the estimate 

g(K) = l + gi K=l-jK, (27) 

with p2 the NNLO RS invariant defined in Eqs.(8). As we shall discuss in 
the next Section the standard a s (/i)BPT is equivalent to the ABPT with 
RS-dependent gi coefficients, and we shall argue that it is likely to give mis- 
leading determinations of A. 

We finally note that although moments of DIS structure functions ap- 
parently involve a factorisation scale and a renormalisation scale, they can 
be directly related to effective charges [7j and ABPT applied [TT]. The n th 
moment of a non-singlet structure function F(x), 

M n {Q) = I x n - 2 F{x)dx , (28) 
Jo 

can be factorised into an operator matrix element and coefficient function 

M n {Q) =< O n {M) > C n (Q, a(fj,),n, M) . (29) 

Here M is an arbitrary factorisation scale, and \x the renormalisation scale. 
The operator O n will have a corresponding anomalous dimension 7^ (a) such 
that 

M dO n ^ ( n) k+1 

n k>0 

where a denotes a(M). The first anomalous dimension coefficient 7q™' ) is 
independent of the factorisation scheme (FS), but the higher coefficients are 
FS-dependent. The coefficient function has the perturbative expansion 

Cn = l + $>^a\ (31) 
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where a denotes a(/i). Integrating Eq.(30) one obtains 

(\ 7o 1° 
exp(J(a))(l + ^rf^). (32) 
' ' ^ fc>i 



A n is an entirely non-perturbative normalisation constant. If one expands 
a(M) as a series in a(/i) one can recast Eq.(32) in the form 



n 7o 7 6 



M n = A n {cR n ]'° , (33) 
where 7?. n is an effective charge, with the perturbative expansion 

n n = a + J2r^a k+1 . (34) 

k>l 

The depend on fi and the non-universal beta-function coefficients, but 
are M- independent, for instance 

r<"> = bin ( J-) - bin ( JL) - * r W + Jjl - a . (35) 

Thus ABPT applies and TZ n = (M n /A n ) bho /c is related to by Eq.(24), 
where r in JC^ S is the NLO coefficient of Eq.(34) with /i = Q. Simul- 
taneous fits for A n and A.jfg have been performed on CCFR data for F 3 
neutrinoproduction DIS moments in Ref . |14j . 



3 Standard RG- improvement in ABPT language 

In this section we wish to clarify the sense in which we claim that ABPT side- 
steps the scheme dependence problem inherent in the standard a s (yu)BPT 
approach. We shall show that the standard approach corresponds precisely 
to ABPT with RS-dependent coefficients. Suppose that we have completed 
a N n LO calculation of TZ, then we can split G(1Z) in Eq.(12) into an exactly 
known piece Q^ n \TZ), containing the known RS-invariant coefficients (i = 
1,2, . . .n — 1), and the unknown remainder 0^(TZ), 

g(n) = g {n \ii) + g {n \n) 

n—l oo 

= (l + 5>**)+£^ fc - (36) 

k=l k—n 
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Standard N n LO a s (p)BPT is exactly equivalent to ABPT with a particular 
RS-dependent choice for the unknown g^ 1 ' coefficients in Eq.(36). These fol- 

— (n) —(ri) 

low on replacing the RS-invariants p n +i, Pn+2, ■ ■ by RS-dependent p n +\, p„+2) 
in which r n+ i,r n+ 2, . . ., have been set to zero in the expressions for the pi. 
Thus at NLO level one has from Eqs.(8) 

-(1) 2 

P2 = —f\C — r 1 

= r 2 c + 4r 3 _ ^ 

Correspondingly, using Eqs.(13) one finds expressions for the gf, 
_(i) (r 2 + r x c) 



9i 



b 

(Arf — ric 2 ) (rf + 2r\c + r\c 2 ^ 



= ^., + vl ...^.. l ., (38) 

At NNLO we would have 

pf ] = 2r\ - 4rir 2 - r\c - 2r lP2 . (39) 

As a concrete illustration let us suppose that r\ = 10, a value which 
is typical of the NLO corrections at the physical scale p = Q for e + e~ jet 
observables such as thrust. Taking Nf = 5 we then find that standard NLO 
a s (p)BPT at the physical scale is precisely equivalent to ABPT with 

Aj^ = Q(0.06415).F(ft) [1 + 29.374ft - 607.92ft 2 + ...]. (40) 

The terms in the square bracket correspond to 1 + (^^(ft). It is seen that 
the coefficients g± are large. In fact at NLO it is straightforward to write 
down the all-orders result for Q^(TZ), in closed form, one finds jH] 



(41) 



If we take ft = 0.05, which is typical for e + e event shape observables 
at Q = M z , then we find that 1 + C? (1) (0.05) = 14.97, to be compared 
with the value G{7Z) = 1 which is all that is exactly known from a NLO 
perturbative calculation. The apparently large higher-order contributions 
to Q{TZ) result from a particular choice of RS, and there is no reason to 
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believe that the unknown coefficients cji really are large. These coefficients 
are simply unknown and are completely independent of the RS. One therefore 
runs the risk with a s (/x)BPT that misleading estimates of A are obtained. 
We believe that a more sensible approach given an exact N n LO calculation, 
is to approximate Q(TZ) by Q^(1Z). If in fact the uncalculated remainder 
Q^ijV) is large then N n LO perturbation theory is not adequate to compute 
the observable. If it is reasonably small then ABPT will give a good estimate 
of the actual Aj^g value. In contrast using the physical scale if r turns out to 
be sizeable is equivalent to assuming values for the unknown higher order 
coefficients which are so large that fixed-order perturbation theory would be 
invalid, for instance with 1Z = 0.05 the first few terms in the expansion of 
G(TZ) in the square bracket in Eq.(40) are [1 + 1.469 — 1.520+. . .]. By studying 
the scatter in the extracted with NLO ABPT from QCD observables one 
is directly learning about the relative size of the unknown NNLO RS invariant 
p 2 . Applying this method to a selection of e + e~ event shape observables in 
Ref. led to the estimate |p 2 | ~ 50. As pointed out in Kef. 15 further 
information about higher-order corrections can also be gleaned from studies 
of the running of observables with energy Q. For instance if p 2 is small then 
the running of the observable should be well-approximated by the NLO result 

<3 H = - wj2(1 + cK >- < 42) 

We cannot resist one further comment. If an observable has been com- 
puted at NLO level so that r\ is known, a part of the next coefficient r 2 can 
be predicted. Rearranging Eq.(8) one finds 

r-i = {r\ - ric + c 2 ) - p 2 • (43) 

The bracketed combination is an "RG-predictable" part of the NNLO coef- 
ficient which is known given a NLO calculation. If the "RG-unpredictable" 
piece given by the NNLO RS invariant p 2 is small and r\ is reasonably large 
then r 2 may be well-estimated by its RG-predictable part. This method has 
been widely applied to make estimates of higher-order coefficients in QCD 
perturbation theory ^H], and these estimates have often proved reasonably 
accurate. We would argue that such estimates are pointless in the sense that 
if r 2 can be accurately estimated in this way then p 2 is small and NLO ABPT 
will potentially give a good estimate of the actual A^j^. ABPT in fact resums 
RG-predictable terms to all-orders. For instance if we make the simplifica- 
tion that c = the RG-predictable estimate for r n is simply r™. Summing 
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this geometric progression to all-orders assuming the one-loop form for the 
coupling a(Q) = l/61n(Q/Ajjy) leads 0111] to the NLO ABPT result 

This all-orders resummation of RG-predictable terms is precisely equiva- 
lent to the complete resummation of ultraviolet logarithms advocated in the 
CORGI approach HQ. 

The first two sections have reviewed and amplified previous work [UJ E] • 
In the remaining sections we turn to some new extensions of ABPT. 



4 ABPT with massive quarks 

If quarks with pole masses {mj} are present then the dimensional analysis 
statement of Eq.(6) becomes 

>m f(K(Q),{ mj }/Q). (45) 



dlnQ 

This is now no longer a separable differential equation, and so it is unclear 
how to integrate it and recover the ABPT relation between TZ and A.jjg. In 
Ref . [2| the generalisation of Eq. (6) to the case of massive quarks has been 
considered. One now has 

m 



dlnQ 



-f <Jl 2 (l + J2fnK n ), (46) 



n>l 



where the f n {{ m j}/Q) ar e RS-invariant coefficients which are Q-dependent 
unlike the pi invariants in Eq.(7) encountered in the massless case. To derive 
expressions for them one can consider the RG equation |T2J 

Here the sum is over quark flavours i = 1,2, ... , Nf. The series expansion of 
the f3(a)&R/da term in powers of TZ must coincide with that in Eq.(7) in the 
massless quark limit when the second term vanishes, in which case — > pi, 

M^T = ~ m2 i l + cU + YPnn n ) . (48) 

n>l 
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With massive quarks the p n have a dependence on {rrij}/Q and are in general 
RS-dependent. The mass derivative in Eq.(47) will have a series expansion 

g^.=h 0i K 2 (l + J2hniR n )- (49) 

1 n>l 

The h n i are again functions of {rrij}/Q and in general depend on the RS. 
They can be easily calculated by simply substituting the perturbation series 
in powers of a for 71 on both sides of Eq.(49), and equating corresponding 
powers of a. One finds [T7j 

dn 

hoi 



<91nm,- 



h Qi h u = g^{r2-rl) (50) 



Substituting these series expansions into Eq.(47), and equating powers of 7Z 
one easily obtains expressions for the f n RS invariants in terms of the p n and 

hoihni |1Z|) 

fo = b + ^hoi 

i 

fofn = bp n + ^ h 0ihni , (51) 
i 

where for n = 1 p\ — c. 

For general renormalisation schemes the h^h ni and p n are separately 
RS-dependent, with the RS-dependence cancelling between the two terms 
in Eqs.(51) to yield RS-invariant expressions for the f n . However, as we 
shall now show, for mass-independent renormalisation schemes such as MS 
where the beta-function coefficients b, c, C2, ... in Eq.(4) are independent of 
the quark masses m,, the h 0i h ni and p n are both RS-invariant. To see this 
consider the equation for the RS-invariant fi. Since b is a universal RS- 
invariant in mass-independent schemes it immediately follows that Yli hoi is 
RS-invariant. It is then clear that for each i individually hoi must be RS- 
invariant since any rrii/p dependence must cancel for each term in the sum 
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separately. Considering the expression for the RS-invariant f\ and noting 
that pi = c is also a universal RS-invariant in mass-independent schemes it 
follows that hoihu is RS-invariant. The crucial point is now that p 2 can be 
related to these h 0i and hoihu RS-invariants, one has 

dp 2 



dlnm, 



= hoihu - ch 0i . (52) 



The fact that the mass derivative of P2 is RS-invariant is sufficient to demon- 
strate that P2 is RS-invariant. Any RS-dependence in p 2 would have to be 
independent of the m 8 , and would hence persist even with massless quarks, 
but this is excluded because p 2 is an RS-invariant in the massless case. Con- 
sidering the expression for the RS-invariant / 2 the fact that p 2 is RS-invariant 
then implies that h Qi h 2 i is RS-invariant. One can then express the mass 
derivative of p 3 in terms of h 0i , h 0i hii, h ih 2 i showing that p 3 is RS-invariant, 
and so continuing in this iterative fashion one establishes that the hoih n i and 
the p n are all RS-invariants in mass-independent RS's. The fact that the p n 
are RS-invariant in mass-independent schemes makes it plausible that the 
ABPT arising from integrating Eq.(45) has the same form as in the massless 
case, with the JC-r, and g,i involving the same combinations of perturbation 
theory and beta-function coefficients. This is indeed the case as we shall now 
show. 

The differential equation in Eq.(45) is no longer separable with quark 
masses and we cannot directly integrate it to obtain the ABPT relation as 
in the massless case. We proceed by noting that from Eq.(9) in the massless 
case we have (assuming MS subtraction) 

hi-rZZ = KkvU + + V K n U n , (53) 

A MS <<- n > 

where the K n are RS-invariant and Q-independent coefficients, the first few 
being 

K=- h ,K- 1 = ± ,K = r - + c\nc , Kl = ^-j,... (54) 

With massive quarks one expects a similar expression but with the Ki now 
having a dependence on {rrij}/Q. To evaluate the Ki with massive quarks 
one differentiates both sides of Eq.(53) to obtain 

i = -f n\i + f 1 n + f 2 'R. 2 + .. .){^- i ^ + k 1 + 2k 2 tz + ...) 
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+ (K']nK+^ + K' + K' 1 K+...) , (55) 

where we have used Eq.(46) and K[ = dKi/dlnQ. Equating coefficients of 
ln7£, TZ^ 1 , 1Z°, 1Z . . ., on both sides of this equation leads , respectively, to 
the following equations 

K' — , K'_\ = , /oK-i +K' = 0, -f K + /o/i^-i + K[ = , . . . . (56) 

In each case one has a first-order differential equation for Ki which needs 
to be integrated f dlnQ to obtain Ki. The boundary condition is that in 
the massless limit (m« — > or Q — > oo) one must reproduce the coefficients 
in Eq.(54). Integrating K' — gives fT = A, where A is an arbitrary 
constant of integration, the boundary condition then fixes A = c/b and the 
massless result K = c/b is reproduced. Similarly K'_ x = integrates up to 
the massless result. Using Eq.(51) for f and setting = 1/6 the next 
equation becomes 

i 

Noting that hoi is RS-invariant and will be a function of the {rrij}/Q we can 
apply the general result that 



E 



ggtW/Q) = dF({ mj }!Q) 
dlnrrii dlnQ 

to recast Eq.(57) in the form 

This integrates trivially to give Kq = (r\/b)+A, where again A is an arbitrary 
integration constant. If we use the notation to denote T\ in the limit of 
massless quarks then we have 

ri=ri 0) +f 1 ({m j }/Q) , (60) 

where f\ contains the {rrij}/Q dependence and vanishes in the massless limit. 
As we discussed above in a mass-independent scheme hoi, the logarithmic 
mass derivative of r±, is RS-invariant, i.e. \i- independent. This implies that r\ 
is also RS-invariant. Thus as in the massless case changes in scale correspond 
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to translating T\ by a constant, and the integration constant A can be fixed 
so that K = (r/b) + cine reproducing the massless result. We now turn to 
the equation for K[. Substituting in the results for /o and fofi from Eqs.(51), 
and setting K = c/b and K_i = 1/b, use of the result in Eq.(58) to change 
from J^iid/drrii) to — d/d\nQ yields 

This integrates up to give K\ = [(r 2 — r\ — r\c) /b} + A, where A is the constant 
of integration. The boundary condition then fixes A = c 2 — c 2 /b to reproduce 
the massless result K\ = (p 2 — c 2 )/b. 

A succinct proof that to all-orders the massless expressions for Ki are 
reproduced with massive quarks is as follows. Instead of differentiating both 
sides of Eq.(53) one can instead differentiate ^ dl ® m . . Equating co- 
efficients on both sides produces a set of conditions equivalent to Eqs.(56) 
with the replacements f f n — > — ^ hoih ni , and where K' n is replaced by 
Si dinm ' us i n S ^ ne result of Eq.(58) one sees that these new conditions 
are precisely equivalent to the original conditions of Eqs.(56) with all terms 
involving the bp n part of f f n removed. This in turn means that this bp n 
subset of terms must separately vanish, but this condition is just the equa- 
tion that determines a given Ki in terms of the p n in the case of massless 
quarks. So, as an example, we can determine K\ by integrating Eq.(61) as 
we have seen. We can also consider the next condition in Eqs.(56) involving 
K' 2 , 

- fofiK + /o/atf-i - /otfi + K' 2 = 0. (62) 

Splitting this equation into the terms involving bpk and h 0i hki and demand- 
ing, as argued above, that each part separately vanishes, gives the two con- 
ditions 

-c 2 + p 2 ~bK 1 = 
y^j-^hQjhu + -h ih 2i - h Qi Ki) + K' 2 = . (63) 

i 

The first condition determines K x to have the same form as the massless 
quark result. In general the bpi subset of terms in the differential equation 
for K' n build the condition satisfied by K n -\ in the massless quark case. No- 
tice that for mass-dependent schemes such as momentum subtraction the 
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above proof fails because the result of Eq.(58) no longer holds since now 
F({mj}/ fi, Q/fi), where /i is the renormalisation scale. Crucially in the mass- 
independent case the RS-invariance of the hoih n i underwrites the validity of 
Eq.(58). 

We turn in the next section to an attempt to resum large infra-red log- 
arithms for e + e~ event shape observables in the ABPT approach, and a 
phenomenological analysis to extract the apparent size of power corrections 
for the measured thrust and heavy jet mass distributions at various energies. 

5 Event Shape Observables in ABPT 

Event shape variables provide some of the most interesting and useful ways 
to confront QCD calculations with experiment (for a recent review see [20 ). 
Their infrared and collinear safety guarantees that they can be calculated 
in QCD perturbation theory, but fixed order calculations describe their dis- 
tributions rather poorly. This situation can be improved by the recognition 
that for a shape variable y, vanishing in the 2-jet limit, large logarithms 
L = \og(l/y) appear at each order of perturbation theory and must be re- 
summed. However, a full understanding of the observed distributions re- 
quires the introduction of large non-perturbative effects, power corrections 
oc e~ h l a ~ A/Q where Q is some hard scale (e.g. the e + e~ centre of mass 
energy). Although the existence of such effects can be motivated by consid- 
ering simple models of hadronization [21] or through a renormalon analysis 
22J, their magnitudes are not at present calculable in a truly systematic 
way from the QCD Lagrangian. Therefore to fit the data we require the 
introduction of either a phenomenological hadronization model or additional 
non-perturbative parameters. Although this hampers attempts to extract re- 
liable measurements of as from the data, it also provides a good opportunity 
to study the IR behaviour of QCD experimentally. For example, in j2H| it was 
proposed to relate the magnitude of the 1/Q power correction to event shape 
means to the average value of a hypothetical infrared-finite coupling. This 
approach leads to predictions for 1/Q power corrections to all event shape 
means in terms of a single additional parameter, a (/•*/)> the zeroth moment 
(i.e. mean) of the coupling at scales < \i < /// (typically fij ~ 2GeV). In 
[2*lj it was shown how this approach could be extended to apply to event 
shape distrubutions. Since then many experimental studies have appeared, 
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fitting event shape means and distributions simultaneuosly for otjj^{Mz) and 
ao(p,i). Generally an approximate (up to corrections ~ 25%) universality of 
the a values is observed, supporting the hypothesis that power corrections 
can be related to a universal coupling in this way. However, all these fits 
use the MS scheme with the "physical" scale choice fi = Q. For the event 
shape means an obvious alternative is to work within the ABPT framework 
described here. This was first carried out for 1-thrust in [23], and some- 
what reduced power corrections were found compared to the physical scale 
approach. This suggests the possibility that the power suppressed effects are 
partly compensating for missing higher order perturbative terms. Recently 
a more extensive analysis was performed by the DELPHI collaboration ^2] > 
taking into account effects arising from the finite bottom quark mass via 
Monte Carlo simulations. They found remarkably small power corrections 
within the ABPT approach, which for many observables were consistent with 
zero. Indeed, the ABPT predictions with no power corrections whatsoever 
gave a better description of the data than the model of |2H1 121] with a uni- 
versal «q- In light of these surprising results, it is interesting to consider 
applying the ABPT method to the event shape distributions. 

In fact, event shape distributions have previously been studied within 
the ABPT framework To do this an effective charge was constructed 
separately for each bin of the data, and NLO QCD calculations were used 
to extract A^r at centre of mass energy Q = Mz- Non-perturbative ef- 
fects were taken into account by using Monte Carlo simulations to correct 
the data back to "parton-level" distributions. This generally improved the 
quality of the prediction, but with this approach it is hard to see whether 
the ABPT distributions prefer smaller hadronization corrections than the 
a(/x)BPT ones. Moreover, even after these corrections were applied there 
were still two kinematical regions where the effective charge ceased to be a 
good description of the data: the 2-jet limit where large logs enhance the 
higher-order perturbative coefficients, and the region (which exists for many 
observables) where the LO result vanishes, causing r\ — > oo. The latter prob- 
lem is hard to address within the effective charge approach, but the former 
problem can in principle be alleviated by introducing a resummation of the 
effective charge beta function. The basic idea is that large logs appearing 
in the perturbation series for the distribution lead to large logs in the per- 
turbation series for p(7Z) JS]. We can then resum this series to arrive at an 
improved description of the data in the 2-jet region of phase space. 

In this section we first show the effect of replacing the hadronization cor- 
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rections of [H| with an analytical power correction ansatz. For simplicity, 
we use a shift in the distribution by an amount C\/Q. This form can be 
motivated by considering simple models of hadronization or through a renor- 
malon analysis [21] and has been found successful phenomenologically (see 
for example J2])- Although better fits are often obtained using the model of 
[23112115 because we are using a different perturbative approximation to stan- 
dard NLO QCD, the subtractions needed to remove double counting will not 
in general be so simple. Next, we outline how we can arrive at an arbitrar- 
ily accurate numerical approximation to the resummed effective charge beta 
function p(7?.)nll by manipulation of the standard NLL resummed expres- 
sions. We then present results showing the effect of using this modified p on 
fits to the distributions of 1 -thrust (r = 1 — T) and heavy jet mass (ph)- We 
use data taken over a wide range of centre-of-mass energies Q = 35 — 189GeV 
(Refs. [26-37]). Lacking information on the correlation between data points 
we have simply combined statistical and systematic errors in quadrature and 
performed a min-x 2 fit, allowing \ 2 to vary by 4 from its minimum to estimate 
a 2a error. This over-simplistic treatment means that our errors cannot be 
considered reliable, however the central values of Aj^g do give an impression 
of the effect of including the power corrections and logarithmic resummation 
into the ABPT framework. 

Before proceeding it is worth commenting on an objection often raised 
to the appliction of the "method of effective charges" to exclusive quantities 
such as event shape distributions. The idea is that the dependence of the 
physical quantity on multiple scales invalidates the dimensional analysis ar- 
gument as presented here in Section 2. However, as pointed out in [TU], this 
is not really the case. Given an observable 1Z = H(Qi, Q2, Q n ) depending 
on n scales, one can simply re-express it as 71 = TZ(Qi, Q2/Q1, Qn/Qi) = 
TZx 2 ,-,x n {Qi)- Here the Xi = Qi/Qi are dimensionless quantities that can be 
thought of as labelling the effective charge which is now a function of one 
single dimensionful scale Q\. We can then write an effective charge beta func- 
tion for this 1Z describing the energy evolution of our physical observable for 
fixed values of the ratios Xj. Of course, this formal manipulation cannot tell 
us whether the p function we arrive at in this way will be well approximated 
by its NLO terms, which is what we require for most phenomenological appli- 
cations, given the current state of the art in perturbative QCD calculations. 
One reason in particular why this might not be the case is if some of the 
Xi become large - typically this leads to powers of large logs Li = log(xj) 
enhancing the coefficients r n in the perturbative series for TZ. These logs can 
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then appear also in the p n and invalidate the simple NLO form for p. To 
solve this we need to apply some form of resummation, which is precisely our 
goal in this analysis. However, as with any perturbative analysis (with an 
unknown remainder function), we can only really determine the accuracy of 
our approximations a posteriori by examining the quality and consistency of 
fits to data. 

Let us now consider the effects of analytical power corrections on the 
results of 0. The procedure used in |S] was to write an effective charge to 
represent the value of the event shape distribution integrated over each bin of 
the data. First, the Monte Carlo program EERAD |3E] was used to compute 
the NLO perturbative coefficients 3 for each bin 



These were then used to write an effective charge, from which a value for 
Ajj-g could be extracted by feeding the data into (J21J). Here, to introduce a 
fit for C\ we drop this "direct extraction" approach and instead perform a 
minimum x 2 fit for Ajjg and C\. For this to work, we need to exclude the 
regions where the EC approach cannot fit the data. For comparison with 
jH], we choose the same ranges selected there (based on the flatness of ri), 
except that the lower end of the range is made proportional to 1/Q when 
looking at data away from Q = Mz- The reason for this is that sub- leading 
non-perturbative effects are expected to become important for y ~ A/Q |2*4] . 
As we are shifting the predictions before comparing to data we require the 
NLO coefficients evaluated for arbitrary bin edges. We have approximated 
these using a set of order 6 polynomial interpolations from the output of 
EVENT2. We have checked by halving the Monte Carlo bin size to 0.005 
that this induces no sizeable error (using the doubled sampling changes the 
best fit values here by less than 2%). 

The results for 1-thrust and heavy jet mass are presented in Fig^ and 
FigO respectively. 

In the case of thrust it appears that the ABPT results prefer larger power 
corrections (and significantly smaller Ajg-g values). For heavy jet mass, the 
situation is similar, although the differences are not quite as extreme. How- 
ever, in both cases we find comparable Ajjg values to those found in |Hj 

3 For our analysis we actually used EVENT2 Jjjj and we have checked that both pro- 
grams give consistent coefficients. 
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Figure 1: 1-thrust: Fits for A.j^g and C\ within the framework of [5| (solid 
ellipse) and standard NLO QCD peturbation theory (dashed ellipse). In the 
latter case the scale is chosen to be fi — Q, and the effect on the central 
value of a change of renormalization scale by a factor of 2 is indicated by the 
arrows. 2a error contours are shown (from allowing \ 2 to vary within 4 of 
its minimum). The fit range is 1 — T = 0.055Mz/Q — 0.23. 
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Figure 2: As Fig[T]but for heavy jet mass. The fit range is ph = 0.035M Z /Q- 
0.2. 
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using hadronization corrected data rather than an analytical power correc- 
tion ansatz. 4 

Now, in preparation for applying our log resummations, we turn to the 
issue of handling exponentiation in the ABPT framework. The typical form 
of an event shape distribution as a double expansion in a and L = log(l/j/) 
is 

= A LL (aL 2 ) + L- x A NLL {aI?) + L~ 2 A NNLL {aL 2 ) + ■ • • , (65) 

a ay 

where the A functions have a perturbative expansion A(x) = Aqx+A\x 2 + . . ., 
and the A are known up to NNLL accuracy. However, some event shapes 
(including 1-thrust and heavy jet mass) exponentiate, meaning we can also 
write 



Ry(y') = 




= C(a7r) exp(Lgi(aiTL) + g2(airL) + ag^aTiL) + 



ag 3 (anL) H ) + D(an, y) , (66) 

where gi and (72 are known 40J. When working with this form of the distribu- 
tion it is conventional to refer to gi as containing the leading logarithms and 
gi as containing the next-to-leading logarithms. C = 1 + 0(a) is independent 
of y, and D contains terms that vanish as y — > 0. These can be calculated to 
NLO by comparison to fixed order results. However, there is no unique way 
of including this fixed order information into the R y (y') (this is the so-called 
matching ambiguity). For example, it is also legitimate to include the C, D 
terms into the exponent (termed "log R matching"), as the difference is of 
order a 3 . Here we choose to define 

Ryiy') = exp(r TZ) , (67) 

so that all the physics is encoded into a single effective charge. This is similar 
to log R matching in that if we re-expand r^IZ in terms of a and L the C and 

4 However, the values of A-p^- quoted in |S] are actually wrongly normalized for two 
reasons. Firstly the factor of (2c/6)< c / b ) ~ 0.85 was omitted, so the results are really 
values for A^-gr. Secondly, the results of EERAD were normalized to the Born cross- 
section (To, whereas the data are normalized to the total cross-section a, and this was not 
taken into account. Mutiplying the EERAD perturbation series by a correction factor 
o~o/o~ = 1 — o.j'K + •• ■ decreases n by exactly 1, increasing the extracted ^-jjg values by 
e 1/b . So the total correction factor to apply to the results of is (2c/b) ( - c/b 'e 1/b ~ 1.11. 
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D functions will clearly appear in the exponent. However, in this approach 
there is no separate matching ambiguity because once we have picked the 
effective charge the inclusion of C and D is automatically determined. In 
particular the full exact NLO coefficient n in a given RS is reproduced if 1Z 
which solves Eq.(16) is expanded in the coupling a for that scheme, thanks 
to the RS invariant which appears on the lefthand side of the equation. 

To perform a resummation of logs in p(7Z) we first note that the r n have 
the form 

r n = rl h L n + rl hh L n ~ l + --- , (68) 
so the structure of the p n as illustrated in Eq. (JHJ) implies that 

Pn = p^L n + pf L L^ 1 + • • • . (69) 

Because L is a logarithm of a physical quantity, this expansion of the p n is 
scheme invariant. We can thus define resummed, scheme invariant approxi- 
mations to p(TZ) as follows 

oo 

Pll (ft) = -bK 2 (l + cK + J2Pn LLn K n ) (70) 

n=2 
oo 

PNLL(^) = -bn 2 (l + C n + J2(Pn LLn + Pn LL L n - 1 )n n ) . (71) 

n=2 

Crucially p^ L is only a function of the r^ L , and p^ LL depends only on the 
r^ L and r^ LL . This means that we can construct pnll given knowledge of 
7?-nll- I n fact, it is a simple matter to arrive at a numerical approximation 
to both pll and Pnll starting from 7?.nll- To do this, we first construct the p 
corresponding to TZnll- Because the NLL parts of this function only depend 
on the NLL parts of TZnll which are exact by definition, p agrees with pnll 
to NLL accuracy. Then, simply truncating p by numerically taking limits 
(L — > oo with LIZ fixed) allows us to extract the LL and NLL terms. In 
some sense the cTZ term is NLL as it is 0(L _1 ) in this limit, but we include it 
in pll as it is obviously present in the full expression, and this avoids having 
to modify (|16j). As an alternative to this numerical procedure we can evaluate 
the first few terms of pnll analytically by series reversion and composition 
(in practice we use this form for TZ < 0.005 to ensure exact cancellation of 
the singularities between the two terms in Q). These converge rapidly to the 
numerical, resummed results. All of our calculations were carried out using 
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the computer algebra system Mathematica jH], allowing the use of arbitrary 
precision arithmetic in taking the L — > oo limit. 

These pnll and pll functions can be used to make predictions for R(y), 
by inserting them in Eq.(16) and numerically solving the transcendental 
equation. It would in principle have been possible to perform correspond- 
ing resummations for the G(7Z) function, and then use the ABPT Eq.(24), 
but construction of p{1Z) is more straightforward. By taking the differ- 
ence in R(y) across the bins in each experimental data set a comparison 
to data can be carried out, including a 1/Q power correction by using 
Rpciy) = Rpriy — Ci/Q). Once again, we need to choose the fit ranges 
so as to exclude the regions where this effective charge description fails. We 
now discuss this choice. 

After exponentiation, the problem as ro — > remains, and in fact for 
thrust worsens; unfortunately this means we need to restrict the fits to 1—T < 
0.18, ph < 0.24 to obtain good fits in the ABPT approach. This suggests 
that writing an effective charge outside of the exponent as in [§] might be 
more sensible; this is certainly possible and leads to a p function where the 
LL, NLL and NNLL terms are known. However this creates problems with 
the resummation. To see this, note that in this case pll consists of just 
the double logarithmic terms containing powers of L 2 R and a cR 3 term and, 
apart from this cR 3 term, is identical to the p induced by the one-loop beta 
function (3(a) = —ba 2 and the double logarithmic distribution 

Ida = d exp f kL 2 a ^ = 2ke- L L(aexp(-kL 2 a)) = 2ke- L L{TZ) , (72) 
a dy dy 

where k is a constant (4/3 for thrust and heavy jet mass). This distribution 
has a peak as a function of a at a max = 1/kL 2 , and so its inverse only exists 
for TZ < 7Z(a max ) = e -1 jkL 2 . As a consequence, for c = 0, Pll{TZ) vanishes 
at this point (where a branch cut starts). Adding the cR 3 term back in, and 
later adding the rest of the NLL terms, does not remove this branch cut. As 
1Z is evolved from Q = oo it increases until it reaches this maximum value, 
and then its evolution becomes undefined. This value turns out to be too 
small to fit the data. One could possibly "switch branches" of p at this point 
and allow 1Z to decrease again, although this would of course still not provide 
a good fit to the data. Note also that this zero of p does not correspond 
to an "infrared freezing" type behaviour because p approaches the zero as 
a fractional power of 1Z — 1Z max - thus the singularity in (|TB|) is integrable 
and the zero is reached after a finite amount of evolution in Q. Accordingly 
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we are forced into using the exponentiated form despite the problems when 
r - 0. 

A lower cut on the shape variables is also necessary. Although we expect 
the behaviour of our distributions to be improved in the 2-jet region thanks 
to the resummation, the onset of non-perturbative effects more complicated 
than a simple 1/Q shift means that we still need to impose a lower cut. 
In fact, even at the purely perturbative level there are problems with our 
resummed results in this region. These are essentially due to the growth of 
1Z as r\ grows large (i.e. as approaches Q). This leads into a breakdown of 
the series of logarithmic corrections to p (as e.g. the NNLL terms are 0(7?.) 
with respect to the NLL ones). In fact, 1Z eventually grows large enough 
that we encounter a branch cut in p which appears due to the branch cut in 
g\. Clearly this behaviour is unphysical, and must be avoided in our fits to 
data. As discussed above we take the lower cut to be proportional to 1/Q 
and obtain good fits for ph, 1 — T > 0.05Mz/Q- Any bins not lying within 
this range have been left out of the fit; a summary of the data we actually 
used is given in Tables 3 and 4. We have also removed the JADE data at 
35 and 44 GeV from the heavy jet mass fits, since its inclusion dramatically 
worsens the fit quality for all the predictions. 

Results of these fits for thrust are shown in Fig|3J and for heavy jet mass 
in FigHJ For comparison, fits using a s (p)BPT to the same data are also 
shown. The fit range for these could in principle be extended as they do not 
suffer from the r — > problem that afflicts the ABPT), however, in order to 
facilitate a direct comparison between the two approaches we have used the 
same fit range for both. The most notable feature of the results is the stability 
of the Ajfg values found within the ABPT framework as we move from NLO 
to LL and then to NLL accuracy, while the fit quality hardly changes. It 
seems that the effect of the leading and next-to-leading logarithms can be 
mimicked by an increase in C\. These Ajg^ values are, however, still smaller 
than the world average. Some examples of the actual distributions are shown 
in FigEland FigJHl 

To investigate the sensitivity of these results to our choice of fit range we 
have redone the fits for a "low" range and a "high" range. The low range was 
determined by decreasing the upper cut until half the bins were excluded, 
and the high range was determined by increasing the lower cut similarly. The 
effects of these changes on the central values of A^r and C\ are shown in 
Table 1 (for 1-thrust) and Table 2 (for heavy jet mass). 

Lastly, we have also considered the so-called "modification of the logs" 
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Prediction 


A/MeV 


d/GeV 


X 2 /dof. 


NLO ABPT 
NLO a s (/x)BPT 


98,116,118 
446,463,517 


1.02,0.85,0.89 
1.47,1.42,1.28 


24/46,59/88,27/44 
26/46,57/88,23/44 


LL ABPT 
LL a s (/x)BPT 


101,119,108 
371,417,478 


0.80,0.63,0.88 
1.12,1.02,0.83 


23/46,63/88,30/44 
24/46,49/88,22/44 


NLL ABPT 

NLL ABPT (mod) 

NLL a s (/i)BPT 


103,107,93 
104,110,95 
200,233,268 


0.50,0.52,0.93 
0.50,0.47,0.90 
0.94,0.79,0.60 


23/46,84/88,34/44 
23/46, 81/88, 34/44 
23/46, 49/88, 23/44 



Table 1. Sensitivity of our fit values for A^r and C\ to the choice of fit 
range for thrust. The fit ranges are 0.05M Z /Q - 0.1 (low), 0.05M Z /Q - 0.18 
(normal), 0.11M Z /Q - 0.18 (high). 



Prediction 


A/MeV 


d/GeV 


X'ldof. 


NLO ABPT 
NLO a s (/x)BPT 


120,114,142 
236,115,221 


1.19,1.22,0.94 
1.65,1.84,1.27 


19/42,50/84,29/41 
19/42,68/84,37/41 


LL ABPT 
LL a s (/x)BPT 


124,123,128 
185,132,146 


1.06,1.07,1.01 
1.39,1.57,1.33 


21/42,51/84,29/41 
19/42,63/84,36/41 


NLL ABPT 

NLL ABPT (mod) 

NLL a 5 (/x)BPT 


125,127,122 
126,128,122 
114,82,67 


0.99,0.97,1.04 
0.98,0.97,1.05 
1.29,1.49,1.62 


21/42,51/84,29/41 
21/42, 51/84, 29/41 
19/42,64/84,37/41 



Table 2. Sensitivity of our fit values for A^- and C\ to the choice of 
fit range for heavy jet mass. The fit ranges are 0.05Mz/Q — 0.12 (low), 
0.05M Z /Q - 0.24 (normal), 0.125M Z /Q - 0.24 (high). 



that is often invoked in studies of event shape variables. This consists of 
modifying L = log(l/y) — > \og((2y max — y)/y) to ensure that the resummed 
parts of the expression vanish at the upper kinematic limit y ma x (which is 0.5 
for both T and p^). The change in central values is shown in Tables 1 and 2. 
One finds that the fitted values change very little. This is to be expected since 
the restricted fit range automatically ensures that the logarithm is essentially 
unchanged in that region. 

6 Discussion and Conclusions 

In this paper we have advocated the idea of ABPT in which QCD observables 
1Z are directly related to the dimensional transmutation parameter of the the- 
ory. Dimensional analysis implies that this relation, given in Eq.(2), holds in 
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Figure 3: Fits to the thrust distribution for and C\. Solid ellipses use 
ABPT, dashed ellipses a s (/i)BPT (with the arrows showing the effect on the 
central value of varying Q/2 < fi < 2Q). The ellipses indicate 2a errors 
generated by allowing x 2 to vary within 4 of its minimum. For a summary 
of the data used see Table 3. 
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Figure 4: As FigJH] but for heavy jet mass. For a summary of the data see 
Table 4. 
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Figure 5: Examples of our best fit NLO ABPT 1-thrust distributions. Red 
is JADE data at 44GeV, green is DELPHI at 91.2GeV and blue is DELPHI 
at 183GeV. 
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Figure 6: Examples of our best fit NLL ABPT 1-thrust distributions. Red 
is JADE data at 44GeV, green is DELPHI at 91.2GeV and blue is DELPHI 
at 183GeV. 
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Experiment 


Q 


Range 


Data Points 


Source 


ALEPH 


91.2 


0.05-0.18 


7 






133 


0.04-0.15 


4 


E3 


DELPHI 


91.2 


0.05-0.18 


10 




21] 




133 


0.04-0.18 


5 




221 




161 


0.04-0.18 


5 








172 


0.04-0.18 


5 




m 




183 


0.03-0.18 


11 


[29 


JADE 


35 


0.14-0.18 


2 


[30, 




44 


0.12-0.18 


3 


[30^ 


L3 


91.2 


0.065-0.175 


4 


ISU 




189 


0.025-0.175 


6 


[32 


OPAL 


161 


0.03-0.15 


6 


eh 




172 


0.03-0.15 


6 


ESI 




183 


0.03-0.15 


6 




m 




189 


0.03-0.15 


6 


[35 


SLD 


91.2 


0.06-0.16 


3 




36 


TASSO 


44 


0.12-0.16 


1 


[37 



Table 3. Summary of the data used in our fits for thrust. 



Experiment 


Q 


Range 


Data Points 


Source 


DELPHI 


91.2 


0.05-0.2 


8 




28 






161 


0.04-0.2 


4 




M 






133 


0.04-0.2 


4 


EH] 




172 


0.04-0.2 


4 






183 


0.03-0.24 


10 




29 




SLD 


91.2 


0.08-0.24 


3 




36 




ALEPH 


91.2 


0.05-0.2 


8 




26 




L3 


91.2 


0.051-0.216 


7 




SI 






189 


0.03-0.24 


14 




32 




OPAL 


91.2 


0.0625-0.2025 


4 




33 






161 


0.0289-0.2025 


5 










172 


0.0289-0.2025 


5 




3S 






183 


0.0289-0.2025 


5 




SS 






189 


0.0289-0.2025 


5 




SS 





Table 4. Summary of the data used in our fits for heavy jet mass. 
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general, but explicit construction of the form of the function $(7£, {rrij}/Q) 
has only proved possible when TZ is simply related to an effective charge. The 
testing of perturbative QCD then reduces to direct extraction of the univer- 
sal parameter Ajfg from the measured experimental data for the observable 
TZ using the fullest information available about the function $. For effective 
charges and with massless quarks we have shown that this function has the 
structure )C^ S TiJfyQilZ), as in Eq.(24). Here the normalisation constant 
JC^( S is obtainable from a NLO perturbative calculation, and is the only part 
to depend on the subtraction convention used to remove UV divergences. 
J-"(TZ) is a universal function given in Eq.(ll), and if we only have pertur- 
bative information our knowledge of Q(TZ) is restricted to the perturbation 
series in Eq.(12), where the coefficients ^ (Eq.(13)) are RS-invariant com- 
binations of a s (yu)BPT perturbative coefficients and beta-function coeffi- 
cients Cj. Given a N n LO calculation the first n — 1 coefficients g\, g 2 , . . . , g n -i 
are known. So at NLO all we know about G{7Z) is G{7Z) = 1. As we showed in 
Section 3 standard RG-improved a s (/i)BPT is exactly equivalent to ABPT 
with a particular RS-dependent choice for the unknown higher-order RS- 
invariants g$ (i > n). With the "physical" scale choice we showed that these 
unknown coefficients can be extremely large, and that therefore the value of 
Ajfg obtained is potentially unreliable. It seemed to us clearly more sensible 
to simply use only the known terms and truncate Eq.(12) after n — 1 terms 
at N n LO. 

We should note that this dimensional analysis motivation was used by 
Grunberg in Ref. [7j to motivate his method of effective charges. Although 
slightly different in motivation the RESIPE approach of Dhar and Gupta in 
Ref.fTU] is exactly equivalent to Grunberg's approach. In the method of ef- 
fective charges one would truncate the effective charge beta-function B(x) in 
the integrand of G(7Z) in Eq.(16) setting B(x) = 1 + cx + p 2 x 2 + . . . + p n lZ n 
given a N n LO calculation. On exponentiating to obtain Q{7V) one would 
then partially include contributions 0(TZ n ) and higher, so the resulting G(7Z) 
would differ by 0(JZ n ) compared to the strict truncation Q^(1Z). The ABPT 
designation is intended to provide a unifying motivation for these various ap- 
proaches, and to reinforce the fundamental importance of the dimensional 
transmutation parameter. 

In Section 4 we considered how Eq.(24) obtained with massless quarks 
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is modified if the quarks have masses. With massive quarks one would nor- 
mally employ Weinberg's form of the RG equation in terms of a running quark 
mass JH] ■ in contrast we have started from the RESIPE formalism results of 
Gupta, Shirkov and Tarasov ^7] in which the physical pole mass is used, these 
are equivalent to the Lie differential RG equation developed by Bogoliubov 
and Shirkov [T^. Our key result is that in mass-independent schemes such 
as MS the pi RS-invariant combinations of pertubative and beta-function 
coefficients are still RS-invariants, which now depend on {rrij}/Q. We were 
able to prove that the ABPT result of Eq. (24) also holds in the massive quark 
case, with {m.,}/Q-dependent coefficients gi, given by exactly the same com- 
binations of pi as in the massless case (Eqs.(13)). 

In Section 5 we extended the direct extraction of A^g^ from e + e~ event 
shape observables of Ref. [§] to include a resummation of large infra-red log- 
arithms L = \og(l/y). One could relate the observable R y (y') to an effective 
charge 71 by exponentiation, R y (y') = exp(r 7t). One could then numerically 
construct Pll{7V) and Pnll(JI) functions by resumming to all-orders the cor- 
responding pieces of p n in Eq.(69). The LL and NLL predictions for R y (y') for 
a given value of Ajjs then follow on inserting these p{7Z) functions in Eq.(16) 
and numerically solving the transcendental equation. To model 1/Q power 
corrections we fitted to a shifted distribution Rpciy) = Rpt(v — Ci/Q)- 
Whilst in principle straightforward a number of complications arose. In 
particular as 1 — T approaches 1/3 the leading coefficient r goes to zero, 
invalidating the effective charge approach. This places a rather stringent up- 
per limit on the fit range. There are also problems in the two-jet region as 

approaches Q, which determine the lower limit. We also noted that one 
cannot directly relate the observables to an unexponentiated effective charge, 
as in Ref. 19], since in that case Pll{TI) has a branch cut, so that 1Z becomes 
undefined. Simultaneous fits for Aj^ and C\ were performed using data for 
thrust and heavy jet mass distributions over a wide range of energies (see 
Tables 3 and 4). The 2a error contours in Ajjg and C\ are shown in Figs. 3 
and 4. NLO, LL and NLL results are shown for both standard a s (yu)BPT 
with physical scale, and for ABPT. For a s (/i)BPT there is a strong decrease 
in Ajjg going from NLO to LL to NLL, whereas for ABPT the fitted A^j^ 
values are remarkably stable. The fitted value of C\ is also somewhat smaller 
for ABPT . We also investigated the stability of the fits to changing the fit 
range in Tables 1 and 2. The ABPT results show significantly more stability 
than a s (p)BPT. 
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We would conclude that, notwithstanding the limited fit range and tech- 
nical complications from which ABPT suffers, there is evidence that the fitted 
power corrections are reduced relative to the standard approach, although 
not as dramatically as in the DELPHI fits of Ref.[T2"] which are consistent 
with zero power corrections. However in that analysis corrections for bottom 
quark mass effects were made, which were not included in our analysis. Event 
shape means have also been measured in DIS at HERA [12] , and it would be 
interesting to perform a ABPT analysis in that well. Unfortunately, 

as we remarked earlier, it is not known how to construct the $ function in 
Eq. (2) for DIS where one has a convolution of parton distributions and hard 
scattering QCD cross sections. One can, however, apply PMS to choose the 
renormalisation and factorisation scales. An analysis along these lines is in 
progress [IS] . 
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